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Abstract 

Using a single-plaquette approximation, novel magnetoinductance effects in Josephson junction 
arrays (JJAs) are predicted, including the appearance of steps in the temperature behavior 
of magnetic susceptibility. The number of steps (as well as their size) is controlled by the 
kinetic inductance of the plaquette whose field dependence is governed by the Abrikosov vortices 
penetrating superconducting regions of the array. The experimental conditions under which the 
predicted effects should manifest themselves in artificially prepared JJAs are discussed. 

PACS: 74.25.Ha; 74.50.+r; 74.80.-g 
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INTRODUCTION 



Many unusual magnetic properties of Josephson junctions (JJs) and their arrays (JJAs) 
continue to attract attention of both theoreticians and experimentahsts (for recent reviews 
on the subject see, e.g. [1-5] and further references therein). In particular, among the 
numerous spectacular phenomena recently discussed and observed in JJAs we would like 
to mention the dynamic temperature rccntrance of AC susceptibility [2] (closely related to 
paramagnetic Meissner effect [3,4]) and avalanche- like magnetic field behavior of magneti- 
zation [5,6] (well described by the theory of self-organized criticahty [7]). It is also worth 
mentioning the recently observed geometric quantization [8] and flux induced oscillations 
of heat capacity [9] in artificially prepared JJAs as well as recently predicted flux driven 
temperature oscillations of thermal expansion coefficient [10] both in JJs and JJAs. At 
the same time, successful adaptation of the so-called two-coil mutual-inductance technique 
to impedance measurements in JJAs provided a high-precision tool for investigation of the 
numerous magnetoinductance (MI) related effects in Josephson networks [11-14]. To give 
just a few recent examples, suffice it to mention the MI measurements [12] on the JJAs of 
periodically repeated Sierpinski gaskets which have clearly demonstrated the appearance of 
fractal and Euclidean regimes for non-integer values of the frustration parameter, and theo- 
retical predictions [13] regarding a field-dependent correction to the sheet inductance of the 
proximity JJA with frozen vortex diffusion. Besides, recently [14] AC magnetoimpedance 
measurements performed on proximity-effect coupled JJA on a dice lattice revealed uncon- 
ventional behaviour resulting from the interplay between the frustration / created by the 
applied magnetic field and the particular geometry of the system. 

By analyzing the influence of MI on magnetization of SIS-type JJAs in the mixed 
Abrikosov state within a single-plaqucttc approximation, in the present paper wc predict 
yet another interesting effect: the appearance of temperature steps in the behavior of sus- 
ceptibility due to vortices induced magnetic field dependence of the kinetic inductance of 
the array. We also discuss the conditions under which the predicted effects can be observed 
experimentally in artificially prepared ordered arrays of unshunted junctions. 
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THE MODEL 



Strictly speaking, to study the influence of MI effects on temperature behavior of mag- 
netization and susceptibility in realistic arrays, one would need to analyze in detail the flux 
dynamics in these arrays. However, assuming a well-defined periodic structure of the array 
(which is actually the case in most experiments [2,6,8,12,15]), to achieve our goal it is suf- 
ficient to study just a single unit cell (plaquette) of the array. (It is worth noting that the 
single-plaquette approximation has proved successful in treating the temperature reentrance 
phenomena of AC susceptibility in ordered 2D-JJA [2,8,15].) Recall that the unit cell is a 
loop containing four identical Josephson junctions. If we apply an external magnetic field 
B normally to the plaquette, then the total magnetic flux $(-B) threading the four-junction 
superconducting loop is given by $(i?) = BS + L(T, B)I{B) where I/(T, B) is the loop mag- 
netoinductance (see below), S the projected area of the loop, and the circulating current in 
the loop reads I{B) — IciT) sin 0(5), where (t){B) = f [n-|- $(i?)/$o] is the gauge-invariant 
superconducting phase difference across the ith junction {n is an integer and by symmetry 
we assume that [8,15] 0i = 02 = 03 = ^4 = 0), Ic{T) is the field-independent critical 
current of the junction, and $o is the magnetic flux quantum. 

In turn, the loop magnetoinductance reads L{T, B) = Lg+Lk{T, B) where Lg = jiQCi is the 
geometric contribution with a being the perimeter of the loop, and Lk{T, B) = fiQX^(T, B)/a 
is the so-called kinetic contribution with A(T, B) being a properly defined London penetra- 
tion depth. In what follows, we shall assume that the array is in the mixed Abrikosov state 
which means that the field dependence of is due to vortices penetrating superconducting 
regions of the array. More specifically, for Bc2 » B > B^, the kinetic magnetoinduc- 
tance follows the hnear field dependence (dictated by the corresponding dependence of the 
penetration depth [16-18]) Lk{T,B) - Lk{T,0) oc B. 

As usual, the net magnetization and susceptibility of the plaquette are given by 
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J{T)[1- cos (j){B)] + 
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is the Hamiltonian of the system which describes the tunnehng (first term) and mag- 
netoinductive (second term) contributions to the total energy of the plaquette. Here, 
J{T) — ($o/27r)/c(T) is the Josephson couphng energy, and V the sample's volume. 



RESULTS AND DISCUSSION 

To properly address the influence of true MI effects on magnetic properties of the array, in 
what follows we assume that for all apphed fields L{T, B)I{B) 3> BS. In this approximation, 
the net susceptibility of the plaquette will depend on applied magnetic field only via the 
universal parameter (5l{T,B) = 27rIc{T)L{T, B)/^o, namely 

X(T,S) c -xo{T) (cos/3^ - 1/3,) (4) 

where Xo{T) = Ic{T)S'^/{n%V), and / = BS/^q is the frustration parameter. 

Let us analyze the obtained results. Figure 1 shows the behavior of the normalized 
susceptibility xl{T.B) = x{T, B) /[xoiT){df'^L/df y] as a function of I3l{T,B) according to 
Eq.(4). Notice the appearance of characteristic minima (steps). They correspond to the 
number of flux quanta that can be screened by the critical currents in single plaquette. 

To further discuss the predicted behavior of susceptibility, we need to specify the explicit 
temperature and field dependencies of /3l(T', B). As usual [8,10], for the exphcit temperature 
dependence of the Josephson critical current 



Ic{T) = /c(0) 



A(T) 
A(0) 



tanh 



A(T) 
2kBT 



(5) 



we will use the analytical approximation of the gap parameter (valid for all temperatures) 
[19], A(r) = A(0)tanh (2.2^^) with A(0) = 1.76kBTc. 

At the same time, as was mentioned before, the temperature and field dependencies of the 
total inductance L(T, B) arc governed by the vortices driven contribution to the kinetic in- 
ductance [16-18], that is Lk{T,B) = Lfc(T, 0)[1 + 5/5ei(T)] with Lfc(T,0) = /xoA^(T,0)/a = 
Lfe(0, 0)/(l-r7r^) and Bci{T) = Bci{0){l-T^ /T^) assuming a two-fluid model expression 
for the temperature dependence of the penetration depth [10]. 

Figure 2 shows the temperature dependence of the normahzed MI induced susceptibility 
x{T, B)/xo{0) for different values of /3l(0, B). Similar to Figure 1, we see the appearance of 
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FIG. 1: The normalized magnetoinductance induced contribution to susceptibility of a single pla- 
quette xl{T,B) as a function of the universal parameter Pl{T,B), according to Eq.(4). 

wcll-dcvclopcd flux-induced temperature steps at some discrete values of the critical param- 
eter /^^(O, -B). The analysis of Eqs. (4) and (5) near Tc reveals that the number of steps 
n{B), their length AT{B) and height Ax{B) depend on the value of /9l(0, B) as follows 



n{B) = -(3L{Q,B), 

TT 



AT{B) 



TT 



and 



Ax{B) = 71 
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FIG. 2: Theoretically predicted dependence of the normalized susceptibility on reduced temper- 
ature for discrete values of the universal parameter Pi^{0,B) (from bottom to top): /3l(0, i?) = 
IT, 27T, Ait, and Gtt. 

These analytical expressions confirm the predicted correlation (seen in Fig. 2) between in- 
crease of the number of steps (as well as decrease of their size) and magnetoinductance 
parameter L(0, B) oc /3l(0, B). 

To test experimentally the predicted here effects, two-dimensional arrays of unshunted 
Nb — AlOx — Nb Josephson junctions with the following working parameters [8,15] can be 
used: lattice spacing a = AQjim (loop area S = a?), critical current /c(4.2i^) ~ 150/xy4 for 
each junction as well as geometric inductance of the plaquette Lg = iiqg ~ GApH, producing 
/?l(0,0) -30. 

In summary, the influence of magnetoinductance effects on the temperature behavior of 
magnetization of the overdamped SIS-type JJAs in the Abrikosov mixed state was studied 
theoretically. Within a single-plaquette approximation, novel phenomenon was predicted 
which should manifest itself through appearance of steps in the temperature behavior of 
susceptibility. The number of steps is totally controlled by the vortices driven field- dependent 
kinetic inductance of the plaquette. 
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